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摘要 : 射电 干涉 仪 通过 获取 同一 天 体 的 信号 到 达 多 台 射 电 天 线 的 时 间 延 迟 来 实现 信号 相关 ， 并 
此 组 成 高 角 分 辩 率 的 射电 望远镜 ， 在 天 体 的 精细 成 图 以 及 高 精度 定位 方面 发 挥 了 重要 作用 。 在 
射电 干涉 仪 观测 中 ， 大 气 中 的 湿性 成 分 一 一 水 汽 ， 造 成 与 天 体位 置 和 结构 无 关 的 额外 延迟 ， 有 具有 
化 快 、 规 律 性 低 的 特性 ， 难 以 建立 精确 的 模型 ， 因 此 需要 对 其 进行 校正 。 水 汽 辐 射 计 是 测量 
气 中 水 汽 含量 的 重要 观测 设备 ， 可 用 于 校正 水 汽 引起 的 时 延 ， 纠 正 射 电 干 涉 的 相位 波动 。 相 比 于 
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其 他 水 汽 探测 方法 ， 水 汽 辐 射 计 具有 更 高 的 时 间 分 辩 率 。 介 绍 了 水 汽 辐 射 计 的 原理 、 研 究 进展 以 
及 在 国内 外 射电 干涉 仪 中 的 应 用 情况 ， 最 后 展望 了 水 汽 辐 射 计 将 来 的 发 展 趋势 。 

X b d: KANT, 射电 干涉 ; KARAKA; 对 流 层 延 迟 
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1 引 8 


射电 干涉 仪 观测 的 电磁 波 穿 过 大 气 层 产 生 时 间 延 人 运 ， 而 对 流 层 中 的 干 性 成 分 和 湿性 成 
分 是 造成 延迟 的 主要 因素 。 其 中 ， 由 于 水 汽 在 时 间 和 空间 上 具有 高 度 可 变性 ， 大 气 中 湿性 成 
pis 分 的 变化 速度 明显 快 于 干 性 成 分 ， 这 使 得 由 水 汽 引 起 的 湿 延 迟 一 直 是 射电 干涉 测量 技术 中 
= 最 严重 的 误差 之 一 ， 对 精细 成 图 和 高 精度 定位 都 将 产生 不 可 忽视 的 影响 ， 在 实际 观测 中 必须 
对 该 延迟 进行 校正 。 

湿 延 迟 与 大 气 中 的 水 汽 含量 直接 相关 ， 水 汽 辐 射 计 (water vapor radiometer, WVR) 是 
通过 测量 大 气 中 的 水 汽 辐射 量 测定 水 汽 及 液态 水 含量 的 一 种 微波 辐射 计 ， 具 有 高 灵敏 度 、 自 
动 化 操作 、 高 时 间 分 辩 率 等 优点 ， 拥 有 非常 好 的 发 展 前 景 。 从 20 世纪 70 FRE, ANEA 
电 天 文 领域 就 开始 使 用 WVR 技术 探测 大 气 中 的 水 汽 含量 " 。 经 过 多 年 的 不 断 改进 ， 水 汽 辐 
射 计 硬件 结构 的 制备 水 平 不 断 提 高 ， 仪 器 的 稳定 性 和 抗 噪 性 愈加 完善 ， 反 演算 法 模型 的 精度 
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也 在 逐步 提高 。 本 文 将 对 近年 来 水 汽 辐射 计 的 研究 工作 进行 简要 介绍 : 第 2 章 介 绍 水 汽 辐射 
计 的 工作 原理 ， 包 括 反 演 大 气 万 线 及 水 汽 te a oe NE 
的 研制 历程 ， 第 4 章 对 水 汽 辐射 计 在 国内 外 射电 干涉 仪 中 的 应 用 情况 作 详 细 介 绍 ， 第 5 章 
结 水 汽 辐 射 计 应 用 于 射电 干涉 测量 时 的 主要 误差 来 源 ， 最 后 为 全 文 的 总 结 及 对 未 来 发 展 方 
向 的 展望 。 


2 水 汽 辐射 计 的 基本 原理 


21 ”获取 天 空 亮 温 度 

大 气 中 的 干 性 成 分 (如 O2, S02, CO, NO, H5S) 和 湿性 成 分 (H20) 都 会 对 电磁 波 产生 吸 
收 作用 。 不 同 物质 对 电磁 波 的 吸收 、 发 射 和 散射 效果 不 同 。 在 微波 波段 ，0。 和 H2O 的 吸 
收 作 用 最 为 明显 。 由 于 水 汽 与 其 他 干 性 气体 成 分 会 造成 射电 信和 号 的 衰减 与 传输 延迟 ， 而 它们 
都 有 各 自 的 吸收 谱 线 ， 通 过 探测 它们 对 应 波段 的 辐射 能 量变 化 ， 就 可 以 对 其 含量 进行 反 演 。 
为 了 便于 表征 和 计算 ， 通 常 使 用 “ 亮 温 度 ” 这 一 物理 量 来 表示 天 线 接收 辐射 能 量 的 大 小 。 物 
体 的 辐射 能 量 用 同 辐射 量 的 黑体 温度 表示 ， 这 个 黑体 温度 值 ， 即 亮 温度 。 射 电信 号 在 穿 过 吸 
收 介质 时 忽略 散射 效应 ， 则 满足 辐射 传输 方程 ， 其 在 天 顶 9 方向 传输 的 解 为 : 


Ta(v) = Tpo(v)e 78°? + seco | kalv, z) T (z)e 7" dz 
> i (1 
T= f ka(v, z)dz 
0 


EH, o AWA, Teo 为 穿 透 介质 前 的 亮 温度 ， 字 为 天 线 接 收 的 亮 温 度 ，ks 为 大 气 吸收 系 
Bl, 2 KAR, T 为 物理 环境 温度 ，r 为 光学 厚度 或 不 透明 度 。 
天 线 接收 到 天 空中 的 微波 辐射 ， 要 先 经 过 信号 放大 器 ， 使 用 转换 器 对 放大 后 的 低频 滤波 
进行 量化 采样 ， 得 到 一 组 采样 电压 ， 即 为 亮 温 度 电 压 。 亮 温度 电压 U 与 天 线 口 面 处 的 目标 
亮 温度 测量 值 Ts 的 关系 可 用 如 下 的 辐射 计 方 程 ”来 表示 : 
U 一 Gs * (Ts + Try) ; (2) 
KH, Gs 为 辐射 计 的 系统 增益 ，TRN 为 系统 等 效 噪声 温度 。 因 此 ， 为 了 保证 亮 温 度 测 量 值 
Ts 的 准确 ， 需 要 先 对 增益 系数 及 系统 噪声 温度 进行 校正 ， 这 个 过 程 称 为 定 标 。 常 用 的 定 标 
方法 有 冷 热 负载 定 标 、 天 顶 定 标 等 。 此 处 以 冷 热 负载 方法 为 例 介绍 定 标 过 程 ， 即 通过 获得 接 
收 机 内 部 的 冷 源 与 热源 辐射 温度 与 输出 电压 ， 建 立 如 下 线性 方程 组 
Unot = Gs : (Thot + Try) 
Ucoia = Gs + (Teoia + Try) 
HF, Trot 和 Twa 为 热 和 冷 源 辐射 温度 ， 为 已 知 值 ; 而 Qio 和 Usoa 对 应 输出 电压 ， 为 测 
量 值 ， 通 过 多 次 测量 ， 即 可 求解 待 标定 的 增益 系数 Gs 与 系统 噪声 温度 TRN， 再 根据 式 (D) 
将 测量 得 到 的 亮 盟 度 转换 为 亮 温度 测量 值 。 


(3) 
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2.2 “水汽 和 液态 水 的 吸收 系数 

水 汽 在 22 GHz 和 183 GHz 附近 均 表 现 有 谐振 带 。 在 干燥 的 条 件 下 ， 使 用 183 GHz 吸收 
线 反 演 水 汽 的 精度 优 于 22 GHz 吸收 线 ”。 但 是 当 大 气 路 径 中 的 水 汽 充足 时 ， 高 频 波段 由 于 
波长 短 ， 水 汽 对 相位 的 影响 容易 过 饱和 ， 因 而 对 183 GHz 频率 的 测定 只 应 用 在 高 海拔 、 极 
地 等 特殊 地 区 。 为 此 ， 本 文 将 重点 介绍 水 汽 吸收 在 K 波段 (20~30 GHz) 的 反 演算 法 。 目 前 
在 该 波段 下 ， 与 数据 拟 合 度 较 好 的 水 汽 吸 收 系数 , 模型 是 van Vleck-Weisskopf (VVW) Bi 
线 模型 ”: 


3.5 
k, =4.5671 x i (22) x e[^140-52)]p y? x 


(à) | 1 | 1 | 
vo / [L(v — vo + Av? (vv + Av? | 
RP, Av 为 吸收 峰 的 半 高 全 宽 ， 定 义 为 : 
300^ %6 300\ 12 
Av = 2.784 x 10 PC) + 4.8p, EJ (5) 


SU, vp = 22.235 10 GHz 为 水 汽 的 谐振 频率 ，w 为 频率 ,TT 为 气温 ，P 为 干 空气 气压 ，p， 
为 水 汽 压 。VVAW 模型 水 汽 吸 收 系数 与 频率 的 关系 曲线 如 图 四 所 示 。 


vvw 曲 线 


A: 277.05 K 
压强 : 1037.6 hPa 
水 汽 压 : 3.228 hPa 


水 汽 吸收 系数 /103km7 


16 18 20 22 24 26 28 30 
频率 /GHz 


ik: 假设 给 定 该 模型 地 面 气象 参数 ， 气 温 为 277.05 K. “EW 1037.6 mbar、 水 汽 压 为 3.228 mbar。 水 汽 吸 
收 系数 的 峰值 v (粉色 虚线 ) 在 22~23 GHz 之 间 ， 接 收 到 的 亮 温 度 在 此 频率 段 最 高 。 


1 VVW 水 汽 吸收 线 


除了 水 汽 的 吸收 系数 曲线 以 外 ， 还 需要 得 到 云 或 圾 中 的 液态 水 的 吸收 特性 。 这 是 由 于 液 
态 水 会 吸收 大 量 的 辐射 ，WVR 接收 到 了 来 自 液态 水 的 亮 温度 ， 但 液态 水 对 折射 指数 的 改变 
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很 小 ， 对 湿 延 迟 造成 的 影响 可 忽略 不 计 。 针 对 阴云 或 有 雾 的 天 气 条 件 ， 为 了 避免 因 液态 水 的 
干扰 而 导致 反 演 结果 产生 偏差 ， 要 根据 液态 水 的 吸收 特性 将 干扰 影响 吻 除 。 液 态 水 的 吸收 系 
数 模型 由 以 下 公式 ”给 出 : 


100-0122(201—T) 
ky =F PE , (6) 
RP, p 为 液态 水 滴 的 密度 ， 和 为 波长 。 
2.3 ”气象 参数 转换 
由 式 (HI) 的 水 汽 吸 收 系数 模型 可 知 ， 吸 收 系数 总 和 环境 气温 TT、 气压 P RKA po 
关 。 当 外 界 环 境 的 水 汽 压 和 气温 已 知 时 ， 可 以 由 理想 气体 定律 来 确定 水 汽 密 度 p, (绝对 
湿度 ): 


_ 217p, 
AC S ) 
则 大 气 层 水 汽 柱 可 降水 量 : 
w=— | pan , (8) 
Pw Jo 


其 中 ，p = 106 g- m? NAAKKA BE. 
若 地 面 气象 的 测量 参数 为 相对 湿度 RR,， 则 可 根据 相对 湿度 的 定义 得 到 Ri 与 p 的 
关系 : 


R, = P" x 100% , (9) 


vs 


其 中 ，pus 为 饱和 水 汽 压 ， 可 使 用 Clausius-Clapeyron 方程 ”来 近似 : 


T ae 25.27 一 273 
ee rc um a 1 
Pus = 6 (5) j d) 


2.4 KORRE 

水 汽 吸收 系数 是 正 算 过 程 的 关键 ， 估 算 水 汽 含量 (precipitable water vapor, PWV) 的 过 
程 就 是 求解 式 (D) 的 亮 温 度 的 大 气 微 波 传输 方程 ， 得 到 水 汽 吸 收 系数 公式 中 的 水 汽 压 等 参 
量 ， 再 由 式 (B) 计算 出 水 汽 含量 的 过 程 。 大 和 气 微波 传输 方程 是 第 一 类 Fredholm 方程 。 由 于 
亮 温度 数据 的 信息 度 不 足 ， 导 致 方程 呈现 病态 ， 不 能 求 出 它 的 解析 解 。 

直观 的 反 演 方法 是 分 层 计 算 ， 通 过 在 不 同 的 高 度 角 下 测量 多 组 亮 温度 ， 得 到 基于 不 同 高 
度 下 的 大 气 参量 方程 组 ， 再 将 各 分 层 的 解 进行 曲线 拟 合 ， 得 到 大 气 参 量 轮廓 线 。 目前 可 使 用 
较为 流行 的 回归 算法 和 神经 网 络 算法 求解 该 模型 方程 组 。 该 方法 的 限制 因素 包括 数据 量 和 
经 济 成 本 。 需 要 用 统计 方法 对 大 量 历史 资料 进行 时 间 空 间 匹 配 ， 调 试 模 型 权重 和 参数 。 大 量 
的 观测 还 需要 使 用 多 个 通道 接收 目标 区 域 辐射 ， 该 观测 过 程 对 技术 要 求 高 ， 成 本 也 较 高 。 为 
此 ， 在 实际 应 用 中 通常 也 可 以 建立 单个 统一 模型 ， 将 式 四) 线性 简化 处 理 ， 由 测量 数据 拟 合 
出 亮 温 度 修正 模型 ， 再 由 修正 的 线性 亮 温 度 大 气 微 波 传输 方程 计算 水 汽 含量 或 湿 延 迟 。 这 种 
建立 简化 模型 的 方法 更 加 经 济 简 便 ， 但 测量 误差 偏差 较 大 ， 以 下 将 重点 介绍 一 些 常见 模型 。 
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2.4.1 双 通 道 主 流 反 演 模型 
水 汽 含 量 与 路 径 延 迟 呈 线性 相关 ( 见 式 (四 ))， 因 此 水 汽 含量 和 路 径 延 迟 均 可 作为 模型 
反 演 目标 。 已 有 研究 给 出 主流 的 双 通道 WVR 线性 模型 . 


Ly = ao + aıTp, + a2Tp, 
e -t-n] 


ag = 

o k a ， (11) 
a v2W, 

_ k 
a2 = iB 


HB, L, ABR HEIR, v, v2 NAS PM, Tei, Te 为 两 个 通道 的 亮 温度 ， 
Ter, Teo 分 别 对 应 两 个 通道 接收 的 经 大 气 衰 减 后 的 宇宙 背景 噪声 温度 ，k = 1.723 x 
107? KK.g-1.m3，W,, To 为 相关 折算 函数 ， 展 开 式 可 见 文献 加 为 了 降低 计算 的 复杂 
度 ， 可 作为 常量 拟 合 处 理 。 

水 汽 含量 和 液态 水 含量 也 可 以 通过 回归 方法 进行 反 演 。Guiraud 等 人 "利用 位 于 美国 丹 
佛 的 探 空 仪 6 年 期 间 的 大 气 数据 (通道 频率 20.6 GHz 与 31.6 GHz 的 WVR)， 根 据 该 模型 的 
思想 给 出 了 如 下 的 系数 修正 值 : 


w = —0.19 + 0.11873, — 0.05607 32, (12) 


| = —0.018 — 0.001147p; + 0.0284Tp2 , (13) 

其 中 ，w 为 水 汽 含量 (单位 为 cm), | 为 液态 水 含量 (单位 为 cm), Tai 和 Tee 分 别 对 应 
20.6 GHz 和 31.6 GHz 通道 的 亮 温 度 。 

在 历史 气象 数据 充分 的 基础 之 上 ， 为 了 提升 反 演 精度 ， 可 加 入 相关 耦合 参数 ， 对 模型 进 
行 改进 。 

(1) 增加 地 面 参数 耦合 

由 于 大 气 剖 面 参 数 可 能 与 地 面 大 气 参数 具有 一 定 的 相关 性 ， 因 此 在 上 述 方法 的 基础 之 
上 ， 可 建立 如 下 多 元 线性 回归 算式 ， 增 加 地 面 坏 境 参 数 带 来 的 影响 : 


w = ag + a1Tp, + deTp, +asPot+aaHot+aslo , (14) 


SEA, Po, Ho, To 分 别 为 地 面 的 压强 、 湿 度 和 温度 ，w 为 视线 方向 的 水 汽 含量 ，w (i—1, 2, 
3, 4, 5) 为 系数 。 

(2) 增加 高 阶 项 耦合 

除了 线性 关系 以 外 ， 大 气 剖 面 参数 与 亮 昌 度 之 间 更 有 可 能 是 复杂 的 非 线性 关系 。 对 此 ， 
HEA ERA EX CD) 的 基础 上 又 增加 每 个 亮 温 度 通道 的 二 次 项 以 及 交叉 项 ， 得 到 非 
线性 回归 模型 : 


w = ao +a1Tp, + a2Tp, +asTe +asTe, +asTe,Tey +asPo +arHo+asTo , (15) 
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JB. a; (i=1, 2, 3,...,8) 为 系数 。 
(3) 引入 大 气 不 透明 度 
一 些 模型 中 使 用 大 气 不 透明 度 r 作为 反 演 变量 ,7 与 亮 温 Te 在 频率 v 处 的 对 应 关 


Tu: = Tg 
=] 0 
rv) E = Tp » 
KP, Tam 为 大 气 平均 辐射 温度 ， 在 一 定 季节 可 视 为 常数 270 K^; Tg, = 2.9 K 为 宇宙 背景 
辐射 。 则 计算 湿 延 迟 模型 可 写成 ; 


Ly = ao +417, + a2T2 , (17) 


KEP, a; (i=1, 2, 3) 为 系数 。 
2.4.2 多 通道 WVR 反 演 模型 

多 通道 式 WVR 反 演 方法 依旧 可 以 采用 双 通 道 模型 的 思想 ， 建 立 线性 相关 关系 。 近 年 来 
扫 频 式 甚至 全 采样 式 WVR. 的 研制 ， 通 过 更 多 更 宽 的 频段 获得 亮 温度 ， 使 得 亮 温度 采样 数据 
更 加 充足 。Tahmoush 和 Rogers” 对 一 种 扫 频 式 WVR 作出 介绍 ， 并 根据 水 汽 与 液态 水 的 
吸收 性 质 ， 建 立 了 如 下 的 亮 温 度 修正 模型 以 提高 分 离 水 汽 与 液态 水 的 精度 : 


Tay(v) =a+bv ch, , (18) 


AF, a,b, c 为 待 拟 合 参数 ，j Nk, 以 22.2 GHz 为 中 心 频 率 的 归 一 化 值 。 系 数 a 用 来 拟 合 
仪器 增益 的 影响 。 由 式 (ER) 表示 的 液态 水 吸收 系数 模型 ， 可 以 看 出 h x v^. Db, EUM 
二 项 bv? 用 以 拟 合 由 液态 水 造成 的 影响 。 
对 于 某 一 高 度 下 的 大 气 薄 层 来 说 ， 考 虑 恒温 恒 压 的 条 件 ， 代 入 式 (IS) 中 的 亮 温 度 修正 
结果 ， 注 层 路 径 湿 延 人 运 d 可 用 下 述 公 式 计 算 : 
3.776 x 105p, Ta, 
EIS ' 
RP, p 为 水 汽 压 ， 其 推导 过 程 可 见 区 可 节 。 
除 此 之 外 ， 一 些 应 用 于 射电 干涉 阵列 的 多 通道 WVR 使 用 了 男 一 种 反 演 思路 。 它 并 未 
致力 于 求解 单个 天 线 的 路 径 湿 延 运 ， 而 是 直接 计算 每 两 个 天 线 间 的 路 径 湿 延 迟 差 异 ， 再 
由 式 (EN) 确定 其 相位 差 。 这 同样 可 以 达到 对 干涉 阵 中 不 同 天 线 之 间 相 位 波动 的 求解 目的 ， 
但 要 求 阵列 的 各 个 天 线 安装 相同 配置 的 WVR。 使 用 该 处 理 方法 的 干涉 阵 包括 VLA (very 
large array), ALMA (the Atacama large millimeter/submillimeter array), ATCA (Australia 
telescope compact array) 等 。 以 ATCA 为 例 进行 说 明 ，ATCA 各 个 天 线 处 安置 了 4 个 通道 
数 的 WVR (filter=1, 2, 3, 入， 某 两 个 天 线 间 的 湿 延 迟 AL, 可 通过 下 式 进 行 计算 ; 


d (19) 


filter=4 AT 
il 
AL, = ` Cw. filter t mer (20) 
filter=1 K filter 


其 中 ， Cw, filter 为 权重 系数 ， K fitter 为 校正 因子 ， 详 见 文献 [L5]. 
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2.5 ”路 径 延 迟 与 相位 校正 
射电 干涉 仪 中 的 单个 天 线 接收 穿 过 大 气 层 的 电磁 波 信 号 如 图 加 所 示 ， 其 治 对 流 层 传播 的 
路 径 差 为 : 


AL-L-G- [ms-G- /mw-Dds+s-G ; (21) 
RP, n 为 折射 率 ，5 为 折射 路 径 ，G 为 不 考虑 折射 效应 的 路 径 ， 工 = f nds 为 信号 传播 的 


光 程 。 定 义 大 气 折射 率 N = 105(n — 1)。 对 于 天 项 方向 的 路 径 而 言 ， 几 何 延迟 项 (S — G) 可 
忽略 不 计 ， 则 可 得 ; 


Aq P = 0 S T Nds . (22) 


地 面 


图 2 ”天线 接 收 穿 过 大 气 层 的 电磁 波 信号 示意 图 
_ 3.776 x 10?p,, 


Aug, A 


WO UI E AAT es BE S DERE n] HEUS A XN, 


T? 
此 ， 式 (Z2) 中 天 项 路 径 湿 项 延迟 项 可 写成 
zeni Dv 
ATzenita — 0.3776 | T5305 - (23) 


对 于 一 层 很 薄 的 大 气 来 说 ， 温 度 与 压强 可 作为 常数 处 理 。 此 时 路 径 湿 延迟 为 : 


Aq = 0.3776 x T .As . (24) 


由 式 (D) 知 ， 衰 减 指 数 项 为 高 阶 项 ， 可 忽略 不 计 。 注 层 亮 温 度 的 辐射 传输 方程 在 天 顶 方 
向 的 取 值 为 : 


AT, =k, T- As . (25) 
XX (2a) 和 式 (四 ) 可 以 得 到 薄 层 处 的 路 径 湿 延 迟 ， 表 示 为 : 


DuA7P 
kole 


Aq = 3,776 x 10" 


(26) 
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若 不 考虑 积分 ， 在 大 气 环境 温度 为 280 K 的 常量 下 ， 由 式 四、 式 (B) 和 式 (Z3) 可 以 得 


到 常用 的 近似 估计 : 
L,m63w . (27) 
由 路 径 延 迟 所 引起 的 相位 变化 为 : 
Ab = 360" (28) 


其 中 ， 和 为 波长 ，Ay (单位 为 度 ) 为 校正 相位 。 


3 水 汽 辐射 计 发 展 简介 


从 20 世纪 40 年 代 起 ，van Vleck ^ ”对 水 汽 和 氧气 在 微波 波段 的 吸收 特性 的 研究 为 大 
气 微波 辐射 理论 商定 了 基础 。1965 年 ，Layson 和 Martin ™ 对 辐射 计 的 微波 折射 校正 理论 
进行 了 系统 的 研究 ， 并 给 出 了 Marcor (microwave atmospheric range correction) 技术 以 反 
演 水 汽 延 迟 量 。20 世纪 70 年 代 ， 美 国 喷气 推进 实验 室 (jet propulsion laboratory, JPL) 证 
明了 水 汽 辐射 计 应 用 于 大 地 测量 领域 的 可 行 性 。 这 种 通过 被 动 式 遥 感 接收 电磁 波 信和 号 来 计 
算 湿 延迟 的 方法 也 开始 应 用 于 射电 天 文 领域 " 。1979 4E, Wa" 提出 用 于 对 流 层 路 径 长 度 校 
正 的 被 动 微波 辐射 计 的 最 佳 频率 对 ， 这 为 双 通 道 设 计 的 水 汽 辐射 计 提供 了 充足 的 理论 文 持 。 
之 后 的 双 通 道 水 汽 辐射 计 的 频率 波段 大 多 依 此 方法 选取 ， 并 被 广泛 用 于 射电 天 文 领域 和 大 
气 科学 等 其 他 领域 。 直 到 目前 ， 依 旧 有 许多 观测 站 在 使 用 双 通 道 式 辐射 计 。 瑞 典 的 昂 萨 拉 
空间 天 文 台 (Onsala space observatory) 设立 了 两 台 双 通道 式 水 汽 辐射 计 Astrid 和 Konrad. 
其 中 Astrid 于 1980 年 5 月 在 哥德堡 - 兰 德 维特 机 场 与 无 线 电 探 空 仪 进行 了 首次 对 比 测量 ， 
Konrad 则 于 2000 年 8 月 在 埃 斯 兰 奇 太空 中 心 (Esrange space center) 的 基 律 纳 (Kiruna) 首 
次 运行 。2017 年 ，Forkman 等 人 ”经 过 测试 表明 ， 即 使 这 两 台 仪 器 已 经 服役 多 年 ， 它 们 在 
4 年 中 的 天 项 湿 延 迟 反 演 结果 与 GPS 数据 的 均 方 根 差分 别 为 0.92 cm, 0.75 cm, FFA EEK IA 
较 好 。 

为 了 进一步 提高 WVR 的 精度 ， 多 通道 式 设计 将 逐渐 取代 传统 的 双 通 道 式 设计 ， 成 
为 目前 以 及 未 来 研制 的 主流 。 美 国 Radiometrics 公司 设计 的 MP 系列 产品 辐射 计 和 德国 
Radiometer Physics 公司 设计 的 RPG 系列 辐射 计 是 目前 使 用 较 广泛 的 两 个 产品 。 以 MP- 
3000 系列 为 代表 的 地 基 WVR 是 国内 应 用 较 多 的 产品 ”， 它 开设 了 35 个 频率 通道 ”， 其 
中 21 个 波段 (22~30 GHz) 通道 ，14 个 V 波段 (51~59 GHz) 通道 ， 检 波 方式 为 串 行 调 
频 ， 成 本 较 低 。RPG 系列 辐射 计 的 最 新 型 号 为 第 五 代 RPG-HATPRO-G5， 采 用 14 个 频率 
通道 ， 其 中 7 个 通道 (51~58 GHz) 用 于 温度 廊 线 的 反 演 ，7 个 通道 (22.24~31.4 GHz) 用 于 
水 汽 廓 线 的 反 演 。 比 较 特别 的 是 ， 它 还 为 极 寒 天 气 地 区 和 高 海拔 地 区 提供 了 183 GHz 波段 
通道 ， 来 提升 低 可 沉降 水 汽 (PWV) 地 区 的 仪器 使 用 性 能 。RPG-HATPRO-G5 采取 了 并 行 
直接 检 波 方式 ， 提 高 了 设备 的 可 靠 性 和 抗 干扰 性 。 图 外 简要 概述 了 水 汽 辐 射 计 的 发 展 历程 。 
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1965 年 ， , JPI 19794, wa MP, RPGA 
Layson 等 将 水 ; 提出 了 双 列 等 多 通道 


提出 Marcor LFA KE 水 汽 辐 tt ae + 
B B 测量 领域 最 佳 Soot DURUM 


3 ”水 汽 辐射 计 发 展 事件 图 


4 水 汽 辐 射 计 在 射电 干涉 仪 中 的 应 用 


下 文 将 从 连 线 干 涉 仪 和 甚 长 基线 干涉 仪 两 个 方面 介绍 WVR 在 射电 干涉 仪 中 的 应 用 情 
况 。 连 线 干 涉 仪 是 指 由 多 元 天 线 系统 组 成 的 干涉 阵列 ， 天 线 间 用 电缆 相互 连接 ， 进 行 相对 较 
短 距离 (WATAK) 的 干涉 测量 。 甚 长 基线 干涉 仪 是 指 利用 几 台 相距 更 远 (如 数 和 干 干 米 ) 的 天 
线 协 同 干涉 测量 的 仪器 ， 仪 器 之 间 不 使 用 电缆 相连 ， 而 是 通过 独立 原子 钟 提 供 频 率 校准 并 
实现 时 间 同 步 ， 各 自 记 录 ， 而 后 进行 相关 处 理 ， 可 以 进行 长 距离 干涉 测量 。 各 干涉 仪 中 的 
WVR .配置 情况 见 表 中 辐射 计 在 进行 测量 任务 之 前 ， 需 要 确定 仪器 测量 的 准确 性 ， 即 对 增 
益 值 的 稳定 性 进行 分 析 ， 并 确定 出 仪器 的 系统 噪声 温度 。 亮 温度 测量 的 准确 性 将 成 为 影响 仪 
器 性 能 主要 因素 ， 评 测 仪器 的 性 能 是 大 部 分 观测 站 的 研究 人 员 必 须 完成 的 重要 工作 。 在 仪器 
性 能 稳定 发 挥 的 基础 之 上 ， 将 对 WVR 的 时 延 改 正 结果 与 其 他 改正 方法 进行 对 比 ， 并 提高 温 

延迟 的 改正 效果 。 


表 1 射电 干涉 仪 中 WVR 的 工作 频段 和 通道 数 


干涉 阵 /天线 ”工作 频段 通道 数 
VLA 22 GHz 5 
ALMA 183 GHz 4 
— 
ERTEK ATCA 22 GHz 4 
SMA 183 GHz 3 
Effelsberg 22 GHz 扫 频 式 
CVN 22 GHz 2/14 
长 基线 干涉 人 
甚 长 基线 干涉 仪 HSN 5 37 a 
= JCMT 183 GHz 3 
= Hs 表 中 工作 频段 指 反 演 22 GHz 和 183 GHz 两 个 水 汽 吸收 峰 附近 的 频段 。 


4.1 ERFA 
411 VLA 

甚大 阵列 (very large array, VLA) 是 位 于 美国 新 墨西哥 州 的 圣 阿 古 斯 丁 平原 上 的 大 型 综 
合 孔 径 射电 望 远 镜 ， 于 1980 年 开始 投入 使 用 。 它 由 28 面 直径 为 25 m 的 射电 望远镜 组 成 。 
其 中 一 面 作为 备用 ， 另 外 27 面 望 远 镜 则 排列 成 Y 字 型 ， 并 可 通过 铁轨 来 移动 单个 望远镜 而 
改变 它们 的 阵型 。 目 前 该 望远镜 阵列 拥有 A, B, C, D 四 种 阵型 ， 不 同 阵型 对 应 的 最 长 和 最 
短 基 线 距离 不 同 ， 来 满足 各 种 科学 工作 的 需要 。 早 期 VLA 天 线 采 用 三 通道 的 WVR 进行 短 
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距离 干涉 的 相位 波动 数据 修正 处 理 ， 该 仪器 在 21.0, 23.5 GHz 处 设置 了 带宽 为 750.00 MHz 
的 通道 ， 在 22.23 GHz 处 设置 了 带宽 为 1000.00 MHz 的 通道 。3 个 亮 温 度 接收 的 天 线 频率 相 
近 ， 为 了 便于 研究 水 汽 波动 与 多 个 亮 温度 的 关系 ，Butler 定义 了 观测 量 AT: 


AT = wT} + WoT + W3T% P (29) 


Hh, w = —0.5, w = 1.0, ws = —0.5 24 — MEME, T1, To, Ts 为 三 个 通道 的 亮 温度 。 
三 个 通道 频率 的 位 置 如 图 团 所 示 。 
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图 4 VLA 三 通道 水 汽 辐 射 计 


2000 年 ，Butler 对 WVR 系统 作出 评测 分 析 ”， 并 简要 讨论 了 该 WVR 系统 在 滤波 
器 稳定 性 、 液 态 水 因素 和 波束 发 散 三 个 方面 可 能 存在 的 问题 。 对 于 50 K 的 系统 温度 而 
言 ，VLA 辐射 计 的 增益 稳定 性 Ag/g 须 好 于 0.001 K 左右。 但 测试 结果 表明 增益 稳定 性 在 
0.002~0.001 K 之 间 ， 未 能 达到 目标 要 求 ， 云 层 中 液态 水 的 干扰 非常 明显 ， 需 要 对 三 通道 的 
测量 数据 作出 合理 的 处 理 ;， 辐射 计 与 天 线 的 波束 发 散 角 约 为 1.4 <?， 这 将 导致 实际 传播 路 径 与 
仪器 测量 接收 的 路 径 不 匹配 ， 产 生路 径 误差 。 

在 2004 年 的 一 次 测试 中 ，VLA R WVR 的 增益 稳定 性 成 功 达 到 了 预期 目标 。 
Chandler “@ A" 对 两 个 VLA 天 线 处 的 三 通道 WVR 进行 测试 ， 并 给 出 这 次 测试 的 结 
果 。 对 于 既定 的 小 于 相位 校正 预期 ， 先 验 计算 得 到 在 系统 温度 为 50~100 K 之 间 ， 系 统 增益 
稳定 性 应 该 为 (2~4)x10-4。 测 试 表明 ， 在 10 s 的 校正 时 间 内 ， 测 试 增益 稳定 性 Ag/g 为 
(4~8)x10 习 ， 达 到 测量 目标 要 求 。 

由 于 更 高 的 分 辨 率 和 灵敏 度 、 更 快速 的 成 像 、 更 广 的 频率 用 盖 范围 需求 ， 工 作 了 20 年 
的 VLA 作出 功能 扩展 ， 升 级 为 EVLA (the expanded very large array). EVLA 的 波段 带宽 
加 大 ， 因 此 支持 了 更 多 通道 系统 的 辐射 计 。Chandler 等 人 ”提出 了 对 原 有 为 VLA 服务 的 
三 通道 WVR 进行 改进 的 项 目 计 划 。 该 计划 打算 设计 并 生产 用 于 接收 EVLA 天 线 开 波段 的 
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集成 型 水 汽 辐射 计 (compact water vapor radiometer, CWVR)。 为 了 提高 仪器 的 反 演 精 度 ， 
将 CWVR 升级 为 五 通道 系统 ， 并 在 内 部 组 件 设 计 上 改 为 集成 型 多 芯片 模块 ， 能 有 效 地 缩小 
WVR 的 内 部 空间 ， 使 CWVR 更 容易 安装 到 天 线 中 。 
2018 年 ，G 记 等 人 ”分 析 并 报告 了 VLA 阵列 CWVR 的 实验 室 测 试 结果 。CWVR 用 单 
片 微波 集成 电路 来 代替 原来 分 散 的 电路 部 件 ， 减 小 自身 的 体积 ， 机 身 更 加 小 巧 ， 并 增加 了 区 
分 液态 水 和 水 汽 的 功能 。VLA 的 波段 频率 范围 为 21.40~24.40 GHz, ifj EVLA 在 波段 
的 频率 范围 为 18.00~26.50 GHz， 波 段 带宽 加 大 ， 因 此 五 通道 系统 的 CWVR ,的 使 用 是 可 行 
的 。 在 仪器 稳定 性 方面 ， 该 报告 表明 其 中 一 个 通道 的 输出 波动 较 大 ， 需 要 维修 校正 。 而 通过 
温度 校正 后 ， 所 有 通道 的 增益 稳定 性 均 小 于 2 x 10-4， 满 足 增 益 稳定 性 的 要 求 ， 可 以 应 用 于 
ngVLA (the next generation very large array). 
4.1.2 ALMA 
地 处 智利 的 ALMA (the Atacama Large Millimeter/submillimeter Array) 阵列 是 多 国 合 
作 的 大 型 天 文 观测 设施 ， 位 于 5000 m 海拔 高度 的 Chajnantor 高 原 上 。 由 于 地 理 条 件 独特 ， 
满足 干燥 特征 ， 因 此 采用 183 GHz 频率 的 水 汽 辐 射 计 来 消除 相位 波动 误差 ( 见 第 区 3 节 )。 
ALMA 辐射 计 拥有 一 套 较 为 完整 的 处 理 系统 ，ALMA 的 相关 研究 人 员 开 发 了 “wvrgcal” 程 
序 ， 可 方便 对 WVR 观测 的 水 汽 数据 进行 离线 批 处 理 ， 并 给 出 了 该 软件 的 设计 以 及 程序 的 
详细 说 明 ”。 基 于 WVR 的 相位 校正 原理 以 及 wvrgcal 中 所 使 用 的 算法 在 ALMA 备忘录 
590 中 给 出 。 
2012 ££, Matsushita 45 AU" fi Fl iz WVR 校正 大 气相 位 ， 并 研究 了 WVR 相位 校正 
在 不 同 频率 、 基 线 长 度 和 天 气 条 件 下 的 有 效 性 。 在 这 次 研究 中 ， 他 们 确认 所 有 基线 的 相 
位 稳定 性 达到 了 ALMA 规范 ， 相 位 稳定 性 通常 能 够 提高 23 倍 ， 有 时 甚至 高 达 7 倍 。 从 
不 同 天 气 条 件 下 的 表现 上 看 ， 当 PWV (precipitable water vapor) 分 别 为 0.5 mm, 0.7 mm, 
1.1 mm, 2.9 mm 时 ， 对 应 不 同 基线 长 度 下 的 总 相位 误差 的 2-p ASD (2 -point Allan standard 
deviation) 值 均 有 所 下 降 ; 但 在 PWV 不 超过 0.3 mm 时 ，WVR 相位 校正 前 后 并 无 明显 差 
异 。 他 们 猜测 ， 在 PWV 过 低 的 干燥 条 件 下 ， 干 性 气体 可 能 成 为 大 气相 位 波动 的 主导 因素 ， 
导致 WVR 的 测量 结果 产生 偏差 。 男 外 ， 即 使 在 良好 的 天 气 条 件 下 或 在 WVR 相位 校正 之 
后 ， 相 位 波动 的 均 方 根 也 会 随 基线 长 度 而 增加 ， 长 基线 的 WVR 校正 效果 不 佳 。 
ri 由 于 干涉 阵列 自身 所 处 高 海拔 地 域 ，WVR 在 干燥 条 件 下 的 性 能 校正 问题 是 ALMA F 
= 涉 阵 列 的 关注 点 。 在 PWV 小 于 2 mm 且 天 空 晴 朗 的 情况 下 ，RMS 相位 应 低 至 20 um 左右 ， 
但 实际 上 只 有 极 少数 达到 此 标准 的 数据 ， 其 残余 误差 项 被 认为 是 由 干燥 的 大 气 环 境 所 致 ”。 
Maud 等 人 下 尝试 在 wvrgcal 程序 中 加 入 经 验 比例 因子 ， 通 过 对 WVR 原始 数据 结果 进行 缩 
放 ， 来 消除 测量 带 来 的 误差 。 利 用 ALMA 射电 干涉 阵 观测 HL Tau 恒星 结构 特征 ， 在 标准 
WVR 改正 和 加 入 比例 因子 后 的 WVR 改正 对 比 见 图 轧 增加 比例 修正 后 ， 环 状 结构 更 加 清 
Wi. MA Bc) 的 流量 剖面 图 可 以 看 到 ， 峰 值 流 量 得 到 增加 ， 修 正 后 局 部 更 加 尖锐 ， 即 亮 环 
和 间 队 的 对 比 度 增加 ， 信 品 比 得 到 提高 。 通 过 对 多 组 数据 总 结 发 现 ，PWYV 较 低 的 数据 (小 
于 1mm) 对 应 着 较 大 的 缩放 比例 值 ， 此 时 采用 比例 因子 的 方法 会 对 数据 的 改进 有 较为 可 观 
的 效果 。 这 表明 该 方法 对 干燥 条 件 下 的 测量 数据 修正 有 效 。 
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轮廓 提取 ， 红 色 为 加 入 WVR 比例 因子 修正 后 的 轮廓 线 ;d) 差异 图 像 。 


图 5 ALMA 观测 HL Tau 恒星 的 原 行星 盘 结构 的 标准 WVR 改正 与 加 入 比例 修正 的 对 比 ™ 


2015 Æ, Hunter 等 人 ”利用 高 精度 点 源 的 观测 得 出 ALMA 阵列 的 几何 天 线 位 置 。 他 
们 使 用 水 汽 辐 射 计 来 校准 每 个 天 线 上 方 的 湿 延 迟 ， 并 对 测 得 的 数据 中 的 水 汽 分 压 和 总 压 进 
行 分 离 。 但 校正 后 的 结果 显示 ， 湿 延迟 改正 后 的 影响 依旧 不 能 很 好 地 抵消 天 线 垂直 方向 上 的 
位 置 偏差 。Hunter 等 人 推测 该 结果 是 由 于 使 用 的 水 汽 反 演 模型 过 于 简单 ， 不 能 反应 出 日 党 


的 湿 延 迟 变化 而 导致 的 。 


2018 Æ, Nikolic 4& AU" 描述 了 整个 ALMA WVR 系统 的 设计 ， 设 计 参 数 的 选择 和 数据 
处 理 策 略 。 他 们 对 WVR 系统 进行 了 初步 测试 ， 结 果 表 明 ， 与 以 前 的 系统 相 比 ，WVR 对 改 
善 相 位 稳定 性 有 了 很 大 的 提升 ， 相 位 校正 系数 ( 即 天 空 亮 温度 变化 与 电子 路 径 变 化 之 间 的 关 
A) 高 达 40 K.mm-l。 测 试 指出 影响 ALMA 处 WVR 精度 的 两 个 因素 一 液态 水 (云雾 ) 和 
与 水 汽 信 号 不 相关 的 残留 相位 误差 。 液 态 水 的 影响 是 WVR 目前 面临 的 关键 难题 。 而 残留 误 
差 的 来 源 尚 未 明确 ， 测 试 表 明 : 外 界 条 件 越 干燥 ,残留 误差 越 明显 ， 基 线 长 度 越 长 ， 残 留 误 


HRK. 
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4.1.3 其 他 连 线 干 涉 仪 

ATCA (Australia telescope compact array) 在 六 个 天 线 处 安置 了 WVR 测量 22.2 GHz 
水 汽 的 辐射 变化 ， 用 于 计算 相位 延迟 。 这 些 WVR 拥有 四 个 通道 数 ， 分 别 设 在 16.5, 18.9, 
22.9 和 25.5 GHz， 带 宽 为 1 GHz。2013 年 ，Indermuehle 等 人 四 对 ATCA 阵列 处 的 WVR 
作出 了 详细 的 介绍 。 通 过 建立 模型 校准 系数 ， 对 亮 温 度 差 进行 加 权 得 到 每 个 天 线 对 之 间 的 
相位 差 。 他 们 将 WVR 对 相位 的 修正 与 插值 法 修正 的 结果 进行 了 对 比 (Wem). HF, e 
(correlation efficiency) 是 插值 残余 相位 与 WVR 残余 相位 的 标准 差 o 的 衡量 值 。 可 以 看 到 ， 
在 编号 为 1 一 10 的 短 基线 中 ， 两 种 方法 相差 不 大 ; 而 在 编号 为 11 一 15 的 长 基线 中 ， 使 用 
WVR 修正 残余 相位 的 方法 优 于 插值 法 。 


表 2 ATCA 每 条 基线 残余 相位 修正 对 比 
插值 修正 WVR 修正 


编号 ”基线 长 度 /m Ae 
C Int €Int OWVR | €WVR 

= 1 92 11.0 0.96 9.4 0.97 0.01 
2 230 10.9 0.96 10.3 0.97 0.00 
3 138 3.4 1.00 9.5 0.97 | —0.02 
4 144 5.1 0.99 13.5 0.95 | —0.05 
5 82 7.5 0.98 11.5 0.96 | —0.02 
6 132 7.3 0.98 11.7 0.96 | —0.03 
7 247 4.1 0.94 16.6 0.92 —0.02 
8 216 5.7 0.99 15.5 0.93  —0.06 
9 240 4.6 0.99 15.3 0.93  —0.06 
10 138 10.6 0.97 13.3 0.95 | —0.02 
11 4 500 47.4 0.50 18.0 0.91 0.40 
12 4 408 41.9 0.59 17.5 0.91 0.32 
13 4 270 41.3 0.60 16.2 0.92 0.33 
14 4 378 45.8 0.53 15.9 0.93 0.40 
15 4 383 39.5 0.62 12.2 0.96 0.33 


SMA (submillimeter array) 使 用 183 GHz 的 三 通道 WVR 对 每 个 天 线 处 的 水 汽 进 行 测 
量 ， 以 用 于 校正 短 距离 干涉 的 相位 波动 。2002 年 对 仪器 进行 检查 ”， 内 部 电子 元 件 均 表现 
稳定 ， 影 响 WVR ,性 能 的 因素 主要 来 源 于 混 频 器 和 放大 器 带 来 的 增益 不 稳定 性 和 接收 器 温 
度 引 起 的 热源 噪声 。Martina 使 用 不 同 的 平滑 算法 尝试 前 弱 这 种 由 仪器 噪声 和 增益 偏差 所 
带 来 的 影响 ”。2015 年 ， 在 由 SMA、APEX (the Atacama pathfinder experiment) 望远镜 、 
SMT (submillimeter telescope) 组 网 的 VLBI WME F, WVR 担任 测量 并 校正 各 站 大 气 
不 透明 度 的 任务 ” 。 
4.2 ” 甚 长 基线 干涉 仪 
目前 的 VLBI 网 尚未 实现 为 所 有 台 站 并 置 水 汽 辐射 计 ，WVR 在 VLBI 观测 中 的 应 用 仅 
在 如 下 一 些 射 电 望 远 镜 站 址 上 做 一 些 研究 。 
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4.2.1 埃菲尔 斯 伯 格 望远镜 

埃菲尔 斯 伯 格 (Effelsberg) 100 m 射电 望远镜 位 于 德国 阿尔 山 (Abr Hills)， 是 世界 上 最 
大 的 全 可 动 射电 望远镜 之 一 。 该 望远镜 安装 了 18 GHz 到 26 GHz 的 扫 频 式 WVR， 随 目标 
源 而 转动 ， 用 于 在 高 频 VLBI 观测 期 间 对 对 流 层 延迟 和 不 透明 度 进 行 校正 。2004 年 ，Roy 
等 人 ™ 详细 介绍 了 该 辐射 计 及 其 用 于 不 透明 度 校 正 的 性 能 。 在 大 气 不 透明 度 的 时 间 序 列 测 
试 数 据 中 ， 将 WVR 的 处 理 结果 与 100 m 射电 望远镜 处 理 结果 进行 对 比 ( 见 图 团 ， 实 线 为 
WVR 在 固定 高 度 角 处 测 得 的 不 透明 度 ， 圆 点 为 Bffelsberg 望远镜 分 别 指向 两 个 不 同 高 度 角 
的 光源 并 利用 特定 程序 计算 得 到 的 不 透明 度 。 图 中 两 种 方法 呈现 出 了 很 好 的 趋势 一 致 性 ， 表 
明 WVR 数据 可 用 于 不 透明 度 校 正 。 


100 m RT RCP 
100 m RT RCP 


大 气 不 透明 度 


13 14 15 16 17 18 19 
UT/h 


6 ”2004 年 2 13 日 大 气 不 透明 度 测量 结果 


2007 Æ, Roy 等 人 ”在 测量 基线 Effelsberg 到 Pico Veleta 间 的 VLBI 相位 时 ， 展 示 
了 一 个 7min 时 间 序 列 的 片段 。 通 过 WVR 的 校正 使 路 径 延 迟 的 RMS 由 1.0 mm 降 到 了 
0.47 mm, Æ 240s 的 尺度 下 的 相关 性 由 0.45 提升 到 了 0.86。 并 将 WVR 与 探 空 仪 、GPS 三 
者 解 算 的 天 顶 湿 延 迟 进 行 对 比分 析 ， 它 们 之 间 的 误差 在 15mm 以 内 。 这 表明 WVR ,结果 可 
于 消除 短期 路 径 延 迟 波 动 ， 改 善 大 地 测量 精度 ， 但 还 不 足以 独立 用 于 VLBI 中 对 流 层 延迟 
的 计算 。 

2012 年 ，Cho 利用 Effelsberg 射电 望远镜 的 WVR 观测 ， 对 VLBI 观测 的 路 径 湿 延 迟 
进行 校正 。 他 结合 欧洲 天 和气 预 报 中 心 的 数值 天 气 模 型 ， 距 Effelsberg 望远镜 较 近 的 5 个 GPS 
站 ， 以 及 探 空 仪 的 数据 ， 与 WVR 校正 的 天 顶 温 延迟 进行 比较 。 对 比 不 同时 段 的 观测 结果 ， 
WVR 的 天 顶 湿 延 迟 与 GPS 和 数值 模型 的 结果 相差 10~50 mm. Cho 将 原始 数据 、WVR 处 
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BASSE. WVR 均值 平滑 后 结果 、WVR. 均值 平滑 加 入 偏 移 修正 结果 等 四 种 不 同方 案 进 行 比 
较 ， 最 后 一 种 方案 的 方差 与 平均 偏差 最 小 。 从 整体 上 来 讲 ， 埃 菲 尔 斯 伯 格 WVR 的 观测 结 
尚 不 能 够 完全 校正 由 VLBI 数据 中 大 气 的 湿性 成 分 引起 的 传播 延迟 ， 这 主要 是 仪器 校准 的 不 
完善 所 致 。 
4.2. "PIS VLBI 网 
中 国 VLBI 网 (Chinese VLBI Network, CVN) 的 各 观测 台 站 均 安置 了 双 通 道 水 汽 辐射 
i. "ede Exi 65 m、 北 京 密云 50 m、 昆 明 凤 凰 山 40 m 和 乌鲁木齐 南山 基地 25 m 射 
望远镜 。 目 的 是 为 了 在 深 空 探 测 等 对 实时 性 要 求 比较 高 的 工程 项 目 中 实现 对 湿 延 迟 的 校 
准 。2018 年 ， 窦 世 标 等 人 ”通过 对 CVN 四 个 并 址 站 GPS 与 WVR 解 求 天 顶 湿 延迟 (ZWD) 
之 间 的 对 比分 析 发 现 ， 昆 明 站 和 乌鲁木齐 站 WVR 安置 区 域 相对 开阔 ， 测 量 数据 系统 差 明 
显 ， 容 易 进 行 GPS 标 校 ， 上 海 站 所 处 地 理 位 置 空气 湿润 ，WVR Bo ZWD 结果 震荡 剧烈 ; 
北京 站 WVR 安置 附近 的 障碍 物 较 多 ，GPS 与 WVR 结果 的 标准 差 较 大 。 上 海 站 与 北京 
的 WVR 数据 系统 差 不 明 显 ， 不 易 进 行 GPS 标 校 。 利 用 2015-2016 年 CVN 四 个 台 站 处 自 
GPS 天 顶 湿 延迟 数据 作为 参考 ， 来 标 校 四 个 并 址 站 WVR 反 演 的 天 项 湿 延 迟 数据 。 经 过 标 
校 后 ， 消 除了 各 台 站 处 的 WVR 数据 与 GPS 数据 的 系统 偏差 (以 上 海 站 为 例 见 图 团 。 除 此 
之 外 ， 他 们 还 对 上 海 站 和 乌鲁木齐 南山 站 的 WVR. 稳定 性 进行 检验 。 上 海 站 WVR 两 个 通道 
本 身 的 亮 温度 振荡 幅度 可 达 2 K。 乌 鲁 木 齐 南山 站 WVR 测量 黑体 的 物理 温度 ， 在 持续 时 间 
10~20 h 处 ， 两 个 通道 的 亮 温度 值 与 黑体 温度 实际 值 偏差 最 大 达到 20 K; 这 表明 仪器 的 观 
测 稳定 性 较 差 ， 需 要 提升 制备 水 平 。 
2020 年 ，CVN 各 天 线 引 入 了 新 的 WVR。 其 中 ， 密 云 站 、 南 山 站 、 昆 明 站 分 别 配 置 了 
新 的 国产 双 通 道 WVR， 改 进 了 天 线 置 的 设计 ， 使 其 不 易 积 水 。 与 原配 备 的 WVR 相 比 ， 亮 
温度 误差 技术 指标 由 1.5 K 提升 为 1 K。 上 海天 马 站 安装 了 一 台 德 国 WVR， 采 用 14 个 通道 
反 演 温度 轮廓 线 和 湿度 轮廓 线 ， 并 提供 四 种 定 标 方法 为 各 个 通道 的 标定 校对 提供 保障 : 黑体 
定 标 、 液 N 定 标 、 噪 声 注入 定 标 、 天 项 扫 角 定 标 。 通 道 数 的 增加 使 测量 效率 提高 ， 在 数据 
量 上 可 以 满足 神经 网 络 算法 确定 回归 模型 系数 的 实现 ， 在 理论 上 可 使 测量 性 能 得 到 显著 提 
升 。 
4.2.3 其 他 天 线 

DSN (deep space network) 天 线 使 用 AWVR (the advance water vapor radiometer) 对 
甚 长 基线 干涉 测量 任务 进行 辅助 校正 。 它 具有 22.2 GHz, 23.8 GHz, 31.4 GHz 共 3 个 频率 通 
道 ，400 MHz 带宽 。2001 年 ，Linfield” 介 绍 了 架设 于 DSN 天 线 副 反射 面 上 的 WVR。 这 是 
一 种 与 天 线 同 轴 转 动 的 WVR， 可 以 有 效 地 减 小 指向 误差 。 文 章 主 要 研究 了 波束 宽度 对 校正 
性 能 的 影响 。 从 结果 上 来 看 ， 其 性 能 由 低 到 高 的 排序 依 此 为 距 天 线 50 m 的 4° FWHM (full 
width at half maximum) 的 WVR， 天 线 同 轴 的 4 各, 2° 和 1° FWHM 的 WVR。 其 中 ， 效 果 
最 佳 的 FWHM 的 WVR 可 以 使 对 流 层 湿 延 退 校 准 的 阿 仓 方才 在 10 s 的 时 间 尺 度 下 降低 至 
校准 前 的 在 100 s 的 时 间 尺 度 下 降低 至 校准 前 的 5p 
JCMT (James Clerk Maxwell telescope) 使 用 183 GHz 水 汽 辐射 计 进行 辅助 校正 。2001 
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天 顶 湿 延迟 /mm 


SHAO 站 2016 年 积 日 /d 


天 项 湿 延 迟 /mm 


YE: a) 标 校 前 ，b) 标 校 后 。 


0 95 
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b) 
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图 7 SHAO 站 GPS 与 WVR 天 项 温 延 迟 (ZWD) 标 校 前 后 的 对 比 


年 ，Wiedner 5& A" 利用 由 JCMT 望远镜 与 CSO (Caltech submillimeter observatory) 所 组 


成 的 干涉 仪 在 0.85 mm 波长 进行 观测 ， 并 通过 两 个 183 GHz 的 WVR 进行 相位 校正 。 结 
果 显 示 ， 通 道 1 和 通道 3 的 相位 校正 效果 较 差 ， 但 通道 2 的 校正 效果 最 优 ， 可 使 相位 波动 


从 141 um 降 到 61 um. 2008 Æ, Dempsey 和 Friberg” 使 用 WVR 对 大 气 校准 的 优化 ， 将 
HARP (heterodyne array receiver program) 计算 出 的 不 透明 度 曲线 与 WVR. 的 不 透明 度 曲 
线 在 高 仰角 处 进行 对 比 ( 见 图 团 ， 结 果 吻 合 度 较 好 ，WVR 可 以 在 短 时 间 的 测量 尺度 上 为 计 
算 大 气 不 透明 度 作 辅 助 优化 ， 提 高 VLBI 任务 中 的 湿 时 延 校正 性 能 。 


5 误差 来 源 


由 于 目前 对 WVR 技术 的 掌握 程度 还 不 算 成 熟 ， 水 汽 辐射 计 反 演 路 径 湿 延迟 的 数据 中 常 


存在 较 大 的 误差 。 造 成 数据 精度 不 足 的 原 


因 可 能 有 以 下 几 个 方面 : 


(1) 液态 水 造成 的 影响 。 对 流 层 的 湿 延 迟 来 自 于 大 气 中 的 水 汽 ， 而 液态 水 对 信号 传播 时 
的 折射 率 影响 很 小 ， 因 此 造成 的 路 径 差 异 很 小 。 但 是 使 用 WVR 接收 天 空中 的 微波 辐射 时 ， 
液态 水 吸收 了 大 量 的 亮 温 度 辐射 ， 并 在 后 续 的 反 演算 法 中 增加 了 其 导致 的 路 径 延迟 。 只 有 将 
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+ Harp skydip 
— WVR | 


大 气 不 透明 度 


20 30 40 50 60 70 80 
高 度 角 /(*) 


ik: 经 过 校准 后 的 结果 在 图 中 用 红色 加 号 表示 ， 实 线 为 同一 时 段 标准 WVR 输出 的 大 气 不 透明 度 结果 。 


图 8 2008 年 5 月 9 日 使 用 HARP skydip 测量 不 同 高 度 角 处 的 大 气 不 透明 度 四 
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测量 得 到 的 水 汽 和 液态 水 的 亮 温 度 完全 分 离 ， 才 能 解决 液态 水 所 带 来 的 影响 。 这 也 是 水 汽 辐 


射 计 反 演 测量 方法 中 的 一 大 难题 。 


(2) 仪器 本 身 的 噪声 。 由 公式 团 可 知 ， 由 WVR 仪器 本 身 的 系统 温度 噪声 TRN 和 增益 


系数 Gs 决定 了 接收 到 亮 温 度 的 准确 性 。 在 比较 干燥 的 条 件 下 ， 亮 温度 值 较 低 ， 上 出 


CL 时 系统 温 


度 噪声 占据 了 主导 地 位 ， 而 增益 系数 的 偏差 会 造成 比例 误差 因子 。 这 都 将 降低 亮 温 度 测 量 值 


的 准确 性 ， 并 导致 最 终 修正 的 路 径 延 迟 偏差 变 大 。 


(3) 大 气 模型 的 精确 性 有 待 提高 。 如 对 流 层 温 湿度 垂直 分 布 模型 存在 误差 ， 液态 水 模型 


不 够 准确 ， 未 考虑 干 性 气体 对 折射 率 产 生 的 影响 等 。 由 不 同 的 温度 和 压强 导致 干 性 气体 对 信 
号 的 折射 率 发 生变 化 ， 从 而 使 路 径 长 度 发 生变 化 引起 的 误差 。Rogers” 提出 这 种 由 气体 中 


的 干燥 成 分 所 造成 的 路 径 变化 可 能 占 到 3096. 


(4) 波束 不 匹配 。 在 相同 的 信号 波长 下 ， 由 于 天 线 本 身 的 口径 不 同 ， 对 应 的 分 辩 率 不 
同 ， 天 线 的 半 功率 波束 宽度 不 同 。 图 日 展示 了 天 线 波束 不 匹配 造成 的 影响 。Tahmoush 和 
Rogers 中 分 析 得 到 ， 对 于 一 个 直径 为 10 m、 接 收 波长 为 3mm 的 射电 天 线 ， 假 设 WVR 的 天 


线 口径 为 60 cm， 那 么 WVR 波束 只 与 该 天 线 在 同 指向 方向 的 400 m 处 完全 匹配 。 


波束 失 配 


WDT WVR 估计 的 湿 延 迟 路 径 变化 与 实际 传播 的 路 径 变 化 的 相关 性 ， 平 均 可 导致 0.3 mm 


的 延迟 误差 。 


(5) 目标 源 的 指向 误差 。 WVR 的 安放 位 置 不 与 天 线 协同 转动 ， 此 时 WVR 的 指向 需要 
与 天 线 指 向 保持 一 致 。 即 使 大 气 较为 均匀 的 条 件 下 ， 指 向 中 的 抖动 也 会 导致 比较 明显 的 路 径 


长 度 变 化 。 在 20° 物理 环境 温度 ，WVR 高度 角 指 向 为 20"，22.2 GHz Ae tba RE d 


FEW 30 KET, 0.1° 的 仰角 误差 会 造成 0.6 mm 的 路 径 长 度 变化 一 。 


C— 
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TE: 橙色 区 域 代 表 只 通过 天 线 采 样 的 部 分 ， 蓝 色 区 域 代表 只 由 WVR 采样 的 部 分 ， 绿 色 区 域 代 表 两 者 都 进行 了 采 
样 。 蓝 色 区 域 的 辐射 信息 是 WVR 的 多 余 信息 ， 橙 色 区 域 的 辐射 信息 是 WVR 的 未 接收 信息 ， 是 由 波束 失 配 引 
起 的 误差 。 可 以 看 到 ， 两 者 只 在 三 种 颜色 交界 的 距离 处 是 完全 匹配 的 。 


图 9 WVR 与 天 线 波束 不 匹配 示意 图 


(6) 反 演 算法 模型 的 误差 。 目 前 使 用 较 多 的 统计 分 析 回归 模型 、 选 代 算 法 模型 等 均 存在 
着 不 同 程度 的 误差 ， 不 能 真实 地 反应 出 亮 温度 与 实际 大 气 轮廓 曲线 之 间 的 关系 。 使 用 更 加 灵 
活 的 神经 网 络 模型 算法 反 演 水 汽 的 垂直 分 布 可 减 小 模型 误差 ”， 但 是 需要 提供 大 量 的 历史 
数据 。 

(7) 其 他 干扰 所 导致 的 误差 。 如 来 自 于 机 箱 顶 部 露水 造成 的 干扰 ， 目 标 源 与 太阳 比较 近 ， 
易 受到 太阳 辐射 的 干扰 ， 目 标 源 处 于 低 仰角 处 ， 容 易 接收 地 面 辐射 等 。 


c 


6 结束语 


本 文 综述 了 WVR 的 仪器 设计 的 发 展 ， 反 演算 法 的 概况 以 及 在 国内 外 射电 干涉 仪 领域 上 
的 应 用 情况 。WVR 的 发 展 已 有 50 多 年 ， 它 在 实时 测量 、 高 度 自动 化 、 高 灵敏 度 以 及 不 受 
地 理 位 置 影响 等 方面 的 便捷 性 和 优势 ， 是 其 他 探测 水 汽 含量 技术 所 不 具备 的 。 目 前 大 部 分 水 
汽 辐射 计 正 在 向 更 多 的 频率 通道 数 发 展 ， 以 得 到 更 加 准确 的 亮 温度 测量 值 。 综 合 来 讲 ， 水 汽 
辐射 计 在 未 来 仍 将 是 消除 对 流 层 湿 延迟 误差 的 研究 热点 ， 具 有 非常 好 的 发 展 潜力 和 前 景 。 

文章 从 原理 上 介绍 了 反 演 路 径 延 迟 的 理论 。 但 在 实际 应 用 中 ， 复 杂 的 气象 条 件 会 导致 诸 
多 误差 因素 ， 多 采用 反 演 模型 计算 。 在 观测 站 放置 WVR 时 ， 需 要 利用 大 量 常规 历史 资料 对 
仪器 和 模型 进行 调试 。 因 此 ， 不 同 的 站 通常 对 应 着 不 同 的 模型 参数 ， 地 域 性 差异 较 大 。 由 于 
精度 具有 偏差 和 不 稳定 性 ， 目 前 WVR 还 不 能 用 于 主动 校正 VLBI 对 流 层 的 湿 延 迟 ， 而 是 多 
以 辅助 改善 湿 延 迟 结果 为 主 。 根 据 近 年 来 人 们 对 WVR 的 研究 ， 还 没有 比较 合适 的 方法 来 解 
决 长 基线 干涉 的 误差 问题 ， 其 发 展 瓶颈 主要 是 液态 水 干扰 项 。 从 硬件 上 可 以 增加 通道 或 扫 
频率 以 提升 拟 合 精度 ， 如 Effelsberg 望远镜 处 的 扫 频 式 WVR, MP 及 RPG 系列 的 WVR 
新 产品 也 拥有 多 个 通道 数 。 从 反 演 函数 的 拟 合 角度 ， 可 以 使 用 神经 网 络 对 可 能 产生 影响 
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参量 自动 归纳 总 结 ， 如 添加 季节 标签 或 年 积 日 参数 ， 添 加 地 面 气象 数据 等 。 在 拥有 大 量 历史 
资料 的 基础 上 ， 可 以 实现 网 络 权重 的 拟 合 。 随 着 人 工 智能 领域 的 迅速 发 展 ， 算 法 的 改进 也 许 
能 在 WVR 的 实 距 上 带 来 更 多 突破 。 另 外 ，Kawaguchi 中 提出 可 通过 对 水 汽 共振 周围 的 频 江 
测量 ， 将 共振 频谱 拟 合 到 理论 值 ， dee d 来 达到 分 离 水 汽 和 液态 水 的 目 


的 。 
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目前 已 有 新 型 数字 水 汽 光 谱 仪 实现 此 方法 ， 然 而 其 在 解决 WVR 的 液态 水 干扰 问题 中 的 


效果 还 有 待 将 来 实测 数据 的 检验 。 
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Application of Water Vapor Radiometer in Radio 


Interferometer 
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(1. Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China; 2. U- 
niversity of Chinese Academy of Sciences, Beijing 100049, China; 3. Shanghai Tech University, Shanghai 
200030, China) 


Abstract: The radio interferometer achieves signal correlation by acquiring the time delays 


of the signal from a common celestial objects between radio telescopes, and thus forms a 


radio telescope with high angular resolution, which plays an important role in fine mapping 


and high precision positioning of celestial objects. However, as a wet component of the at- 


mosphere, water vapor in radio interferometer observation results in extra delay independent 
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of the position and structure of celestial bodies, which has the characteristics of fast varia- 
tion and low regularity, hence it is difficult to establish an accurate model. Accordingly, the 
delay due to water vapor needs to be corrected. Water vapour radiometer is an important 
observation equipment for measuring water vapour content in the atmosphere. It can be 
used to correct the time delay caused by water vapour to reduce the phase fluctuation of 
radio interference. Compared with other water vapor detection methods, the water vapor 
radiometer has higher time resolution. In this paper, the principle, research progress and 
application of water vapour radiometer in radio interferometers are introduced. Finally, the 


trend of water vapour radiometer development in the future is prospected. 


Key words: water vapor radiometer; radio-interferometer; atmospheric profile inversion; 


troposphere delay 


